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The effect of random shooting of particles is considered on the basis of so- 
lution of the Schrodinger equation and in terms of the Wigner function. Two- 
particle description shows, in particular, that initial correlation leads to high 
velocities of particles. This could be a potential mechanism for obtaining energy. 
Evolution of the n-particle probability distribution is described analytically. 

The Heisenberg's uncertainty relation reads: 

ApAx ~ h (1) 

where Ap is uncertainty of momentum {p — mv - product of mass and velocity of 
particle). Ax is uncertainty of coordinate and h is Planck's constant. Relation 
(1) seems to preclude the existence of trajectory: the more exactly we know 
location the more uncertain is velocity. Probability density of particle position 
is given by P = \[/\l>* ^ where ^^(tja:) is the wave function and star indicates 
complex conjugate. For simplicity, we start with the one-dimensional case, but 
the 3-D generalization is simple (see below). 

Schrodinger's equation has the form (see, for example, Ref. 1): 

iW'^/dt ^ -(hV2m)a2*/aj;2 + [/(a;)^' (2) 

Here U{x) is the potential energy of the particle in the external field, i 
is the imaginary unit. The strange behavior of particles is not associated with 
particular form of potential energy. So, we put U = and consider the influence 
of the quantum vacuum on the motion of particle. Equation (2) becomes very 
simple: 

d'i'/dt = iqdH/dx^, q = h/2m (3) 

where q is the coefficient of imaginary diffusion (ID). Equation for has the 
same form as (3) with minus in the right hand side. 

General solution of equation (3) in unbounded space can be obtain with use 
of Fourier transform: 

d^/dt^ -iqk"^^ (4) 

where k is the wave number and 5'(<, k) is the transform of ^'(i, x). Solution of 
(4) is: 



*(t, k) = ^a{k) cicpi-iqk^t} 



(5) 



where zero indicates initial function. The inverse transform of (5) gives: 

^{t,x) = {'imqt)-^'^ j d:c'*o(a;')exp{-^^^^^} (6) 

— oo 

Wc would like to clarify the physical sense of ID and expressions (5) and (6). 
Let us introduce Gaussian random velocity v{t) with zero mean (< v{t) >= 
0) and correlation proportional to 5-function (see, for example, textbook [2]): 

< v{t + T)v{t) >= qd{T) (7) 

Here and below brackets <> indicate statistical averaging. Consider trajec- 
tory, produced by such velocity: 

t 

z{t) =zo + j v{t')dt' (8) 



We assume that initial position zq has Gaussian distribution, independent 

of f (t), with zero mean (< zq >= 0) and dispersion < (zq)^ >= o?- Then, z{t) 
will have Gaussian distribution with zero mean and dispersion: 

< W)? >= + <fi (9) 

Here we used (7) and independence of from zq. 
Now, consider complex trajectory: 

t 

y{t)=yo + {l+i) J v{t')dt' (10) 



Here yo is real with the same statistics as zq. We have < y{t) >= and 

< [y{t)]^ >^ + 2iqt (11) 

High order statistical moments of y{t) can be calculated from these first two 
by the same formulas as for the Gaussian distribution. Effective probability 
density (epd) for yi = y— < y > 

{2n<yl »-V2exp{-^^^} 

can be used, as if y is real. For statistical moments of y*{t) in all formulas 
we have to change the sign in front of i. We will call such complex random 
processes complex-Gaussian. 

Returning to (6), for simplicity of calculations, let us choose: 

^o{x) = (2™2)-i/4exp{-a;V4a2} (12) 
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This wave function corresponds to the Gaussian probabihty density for the 
initial position of particle: 



Po{x) = (27ra2)-V2 ex.p{-x'^ /2a'^) (13) 
Substitution of (12) into (6), after simple calculation, gives: 

= (2.)V4(l2+,^,)V2 ^>^Pi-4(^2^} (14) 

This expression (apart from normalizing factor) is the epd for the complex- 
Gaussian process like (10) with a replaced by v^a. It suggests that quantum 
particles have hidden complex trajectories. It also suggests that the quantum 
vacuum interacts with particles by producing tachyonic impulses (with imagi- 
nary energy), which leads to ID (see below). 

Using (14), we calculate the probability density for the position of particle 
in real world: 

1 a;2 

[27r(a^ + w t )J ' 2(a^ + ■w'^t'^j 

Here w = qa~^ is effective constant velocity. The trajectory for (15) is 
random shooting: 

x{t) = Xq +ut (16) 

where xo and u are independent random constants with Gaussian distribu- 
tions and < Xo >= 0, < (xq)^ >= a^, < u >= 0, < u'^ >= vr^. The fact that 
w is inversely proportional to a is, of course, manifestation of the uncertainty 
relation (1). 

Note that if we multiply the initial wave function (12) by exp{i7a;}, where 
7 is real constant, the initial probability density (13) will not change. However 
in (15) X will be replaced by x — 2q'yt, so we will have < x >= 2q'yt. In this 
paper we choose initial conditions such that the mean positions of particles will 
not change in time. 

The quadratic time-dependence of dispersion in (15) shows that quantum 
complex trajectories like (10) can not be just projected on real axis, leading to 
linear time-dependence (9). The nonlinear operation in (15) is proportional 
to the joint epd for the complex and complex conjugate trajectories, as if they are 
independent. This indicates an interesting sort of interaction. It is remarkable 
that such interaction produces random shooting (16). Note that the effective 
energy mw^/2 is inversely proportional to m. 

Let us look into details of this interaction by considering the characteristic 
function: 



$(i, k) = j dxP{t,x)ex.p{ikx} = j dx'^{t,x)^*{t,x)ex.]){ikx) 



(17) 
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Using Fourier transformation, we get from (17): 



^{t^k)^— / dm'^{t,m)'^*{t,m-k) 




(18) 



The transform of (12) is: 



*o(fc) = 23/47rV4aV2exp{-a2fc2} 



(19) 



Solution (5) gives: 



(20) 



From (18) we now have: 



$(t,A:) = 2i/27r-V2a / dmexp{-a'^[m'^ + {m- kf]- iqt[m^ - {m- kf]} (21) 



which corresponds to the transform of (13). 

Let us stress the difference between the ordinary diffusion, when iq = \ is 
real, and our case of ID. For real A in the second square bracket in (21) wc will 
have sign plus between and (m — fc)^ - the same as in the first square bracket. 
For real A we will have solution of the original equation (3) for the probability 
density (instead of wave function or cpd) with the linear timc-dcpcndcncc of the 
dispersion. In our case of ID the sign inside the second square bracket is minus. 
This leads to elimination of the linear time-dependence and to the solution: 



which is the transform of (15). The elimination effect can be considered as 
interaction between m-tachyon and (m — fc)-antitachyon. This interaction we 
sec after averaging. It will be very interesting to investigate in future - what 
kind of " battle of tachyons" takes place before averaging. 

There is another way of looking at the random shooting of particle with the 
characteristic velocity w = qa^^. Consider a pair of vortices with intensities q 
and —q, placed on the line perpendicular to real line in the complex plane: one 
above the real line, another below on the same distance ^ a . These vortices 
will move each other and the particle with velocity ^ w. This picture can be 
considered as coherent structure, produced by the flow of tachyons. 

The 3-D generalization of ID is natural to do in terms of characteristic 
function with normalization condition <I>(t, 0) — 1. For simplicity, we assume 
isotropy of initial probability. Substitution of 3-D Laplacian A = /dx\ 
-\-d'^ / dx2+d'^ / dx\ instead of d"^ jdx^ in equation (3) will not change the isotropy. 



For t = : 



$(0,A:) = cxp{-- 




(22) 



$(t, k) = exp{--(a2 4- w;2t2)fc2| 



(23) 
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It means that in expressions (12) - (15) wc have now x? = xf+x'^+x'^ and prefac- 
tors are rased to the third power. In equation (23) we now have = kl+k^+k"^- 
Correspondingly: < Xb >= 0,(6 = 1,2,3), < XbXc >= {a^ + w'^t'^)6bc, where 
5bc is the unit tensor. So, wc get 3-D Gaussian distribution for the position of 
particle with indicated moments. 

Tachyonic impulses can be at the core of the phenomena of quantum entan- 
glement, as indicated before [3]. Consider system of n particles, positions of 
which X = {xi,...,Xn) are initially correlated. We start again with 1-D case and 
3-D generalization will follow. We assume that particles do not interact directly. 
The first possibility is that particles may not interact because of their nature. 
The second possibility is that we separate them initially, so that probability 
of direct collision is low (here 3-D generalization is handy). We choose initial 
condition such that mean positions < x > will not change in time (see below). 
To simplify notation, we assume that < x > is already extracted from x. For 
the corresponding wave function ^{t,x.), instead of (3), we now have equation: 

a^=^E«ia^'«i-^ (24) 

where summation is over j = 1, n and mj are masses of particles. Fourier 
transformation of (24) gives: 

d^/dt=-i[J2Qjk]]^ (25) 

Here ^'(i, k) is the transform of ^'(t, x) and k = {ki,...,kn) is the wave 
number vector. Solution of (25) is: 

$(i,k) = $o(k)exp{-zt[^g,fc2]^ (26) 
We assume initial Gaussian distribution for positions of particles [2] : 

^^^'^^ = (2.)V.(ktC)V. exp{-lxC-x} (27) 

Here C =Cji =< XjXi > is the initial covariance matrix (symmetric) and 
C~^is the inverse matrix. In particular, Cn =< xf >= af is dispersion for 
position of the first particle with zero mean (< xi >= 0), aj corresponds to the 
particle j, Pji{0) = Cji{ajai)~^ (with j ^ I) is the initial correlation coefficient. 
Characteristic function for (27) is: 

$o(k) =< exp{ikx} >= exp{-ikCk} (28) 

1/2 

Choosing initial wave function ^'o(x) = [-Po(x)] , after simple calculation, 
we get transform: 

$o(k) = 23"/V"/^(£)eiC)i/^ exp{-kCk} (29) 
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Formulas (26) and (29) give: 



k) = 23"/V"/^(£>eiC)i/* exp{-k[C + itql]k} (30) 

Here ql stands for diagonal matrix with elements qj . This expression (apart 
from constant prefactor) is characteristic function for the joint epd of n complex- 
Gaussian trajectories like (10) with statistically independent velocities Vj{t). It 
means, in particular, that < vi{t + T)v2{t) >=< vi{t + r) >< V2{t) >= and 
mutual correlation function of these trajectories is equal to correlation of initial 
positions. Such simple is the quantum entanglement in the "complex world". 
In the real world situation is more complicated, as we will see below. 

Consider characteristic function for P(t, x) by using generalization of (18): 

*(i, k) = J dm^it, m)$*(i, m - k) (31) 

Substitution of (30) into (31) gives: 

$(t,k) = {-)'^/^{DetCy/^ j dmexp{-2mCm+2mCk-2iimx+zikx-kCk} 

where we used symmetry of C and introduced vector x = (gifci, ...,qnkn). By 
substitution mx = mCC~^x, the integral takes standard Gaussian form: 

i>(t,k) = (-)"/2(£)efC)^/2 / dA*exp{-2/xC/x+ ^xCx + iikx-kCk} 

TT J 2 

where /x = m — x/2, % = k — itC~^>r. The linear timc-dcpcndence of 
dispersions is eliminated by the interaction of the tachyonic impulses, as can be 
anticipated. But, there are interesting details, which are hard to predict without 
calculations. The integration gives characteristic function, which corresponds 
to Gaussian distribution: 

$(t,k) =exp{-^[kCk + iVC-V]} (32) 
For n= 1, formula (32) reproduces (23). For n = 2 formula (32) gives: 

$(t,k) = exp{-^[al{t)kl + 2p{t)ai{t)a2{t)hk2 + al{t)kl]} (33) 



al{t)=al+wlt'^, al{t)=al+wlf, Wi = qi{aiPo)-\ W2 = q2{a2po)~'^ (34) 

... aia2 - wiW2t'^ 
Pit) = Po ^ f,. (35) 

Here uf and are corresponding dispersions, p is correlation coefficient 
with initial value pg /^o = ~ PoY^^- The first important feature of 
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this distribution is the factor /3q in (34) - increase of velocities due to initial 
correlation. In terms of vortices (see above), it could mean such interaction 
that, in average, gets them closer to real line. 

Additionally, wc get evolution of correlation coefficient p{t) from initial value 
pQ to asymptotic value — Pq (when t ^ maxjai/wi, 02/^2}). We can say that 
quantum vacuum "does not like" the measurement- imposed correlation and 
makes the inversion. 

It will be interesting to play with these features experimentally. Particu- 
larly, it is important to investigate if the imposing of initial correlation can be 
made energy-efficient. Then the shooting effect with high velocities could be a 
potential mechanism for obtaining energy. 

The 3-D generalization for the two-particle system is as follows. Positions of 
particles and wave numbers wc denote by .Tj ^and fcj.5, where j = (1, 2) indicate 
particle and b = (1, 2, 3) corresponds to 3-D vector. Matrix of initial correlations 
now has two sets of indexes: Cji^bc- In simplest case: Cjific = CjiSbc- In this 
case we can use (33) with products kf, fcifc2 and fc| replaced by the scalar 
products of vectors, in particular, kik2 is replaced by fci, 1^2, 1+^:1,2^:2, 2 + ^1, 3^2,3- 
From such generalization of (33) it follows: < Xj^i, >= 0, < Xi^xi^c >= ^i^bc 
, < a;2,b.T2,c >= o-lSbc, < a;i,b.T2,c >= /0CTicr2^fec- 

In general case formula (32) gives: 



Velocity of the shooting in (37) is determined by the term of the inverse 

covariance matrix, which is multiplied by ^ in the initial probability density 
(27). For the two-particle system this term, according to (34), is proportional to 
(1 — Po)~^. For n > 3 velocity of the shooting depends of all initial correlation 
between particles. 

The presented results can be also obtained and interpreted in terms of the 
Wigner function (see, for example, recent book [4] and references therein): 



Here we consider the 1-D case for one particle and p is additional (momen- 
tum) variable. The probability density is given by integral: 



< XjfiXl^c >= Cjl^bc + t'^QjqiGji],c 



(36) 



In particular: 



< a;?_i >= Cii,ii +t^giCii^ii 



(37) 





(39) 



Solution (5) [or (6)] corresponds to: 



W{t,x,p) = Wo{x-^,p) 



(40) 
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Initial condition (12) gives: 



^o(x,p) = -exp{---^} (41) 
Using (40), we have: 

Wit..,,) -^exp{- ^ ^ (42) 

Integration (39) of (42) reproduces (15). In addition, we get interpretation 
of the random shooting (16) in terms of variable p. Let us note that formula 
(15) for one particle was known before (see references in [4]), but interpretation 
was different. Expression [a? + w^t^Y/'^ (in our notation) was interpreted as 
trajectory [4]. In our interpretation, the motion is with constant (random) 
velocity (16). 

The ri-particle system has analogous description by the Wigner function 
with interpretation of general formulas (32) and (36) in terms of corresponding 
p- variables. It is helpful to have in mind both approaches ( Schrodinger and 
Wigner ) in order to develop intuition, which unites tachyonic impulses and 
p- variables. 

We plan detailed investigation of systems with three or more particles by 
using presented above general results. 

Further studies of ID can also open new perspective in quantum field theory. 
I thank V. I. Tatarskii for useful discussion. 
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